home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / chegst.z / chegst
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCHHHHEEEEGGGGSSSSTTTT((((3333FFFF))))                                                          CCCCHHHHEEEEGGGGSSSSTTTT((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CHEGST - reduce a complex Hermitian-definite generalized eigenproblem to
  10.      standard form
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CHEGST( ITYPE, UPLO, N, A, LDA, B, LDB, INFO )
  14.  
  15.          CHARACTER      UPLO
  16.  
  17.          INTEGER        INFO, ITYPE, LDA, LDB, N
  18.  
  19.          COMPLEX        A( LDA, * ), B( LDB, * )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      CHEGST reduces a complex Hermitian-definite generalized eigenproblem to
  23.      standard form.
  24.  
  25.      If ITYPE = 1, the problem is A*x = lambda*B*x,
  26.      and A is overwritten by inv(U**H)*A*inv(U) or inv(L)*A*inv(L**H)
  27.  
  28.      If ITYPE = 2 or 3, the problem is A*B*x = lambda*x or
  29.      B*A*x = lambda*x, and A is overwritten by U*A*U**H or L**H*A*L.
  30.  
  31.      B must have been previously factorized as U**H*U or L*L**H by CPOTRF.
  32.  
  33.  
  34. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  35.      ITYPE   (input) INTEGER
  36.              = 1: compute inv(U**H)*A*inv(U) or inv(L)*A*inv(L**H);
  37.              = 2 or 3: compute U*A*U**H or L**H*A*L.
  38.  
  39.      UPLO    (input) CHARACTER
  40.              = 'U':  Upper triangle of A is stored and B is factored as
  41.              U**H*U; = 'L':  Lower triangle of A is stored and B is factored
  42.              as L*L**H.
  43.  
  44.      N       (input) INTEGER
  45.              The order of the matrices A and B.  N >= 0.
  46.  
  47.      A       (input/output) COMPLEX array, dimension (LDA,N)
  48.              On entry, the Hermitian matrix A.  If UPLO = 'U', the leading N-
  49.              by-N upper triangular part of A contains the upper triangular
  50.              part of the matrix A, and the strictly lower triangular part of A
  51.              is not referenced.  If UPLO = 'L', the leading N-by-N lower
  52.              triangular part of A contains the lower triangular part of the
  53.              matrix A, and the strictly upper triangular part of A is not
  54.              referenced.
  55.  
  56.              On exit, if INFO = 0, the transformed matrix, stored in the same
  57.              format as A.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCHHHHEEEEGGGGSSSSTTTT((((3333FFFF))))                                                          CCCCHHHHEEEEGGGGSSSSTTTT((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      LDA     (input) INTEGER
  75.              The leading dimension of the array A.  LDA >= max(1,N).
  76.  
  77.      B       (input) COMPLEX array, dimension (LDB,N)
  78.              The triangular factor from the Cholesky factorization of B, as
  79.              returned by CPOTRF.
  80.  
  81.      LDB     (input) INTEGER
  82.              The leading dimension of the array B.  LDB >= max(1,N).
  83.  
  84.      INFO    (output) INTEGER
  85.              = 0:  successful exit
  86.              < 0:  if INFO = -i, the i-th argument had an illegal value
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.